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Abstract 

The nature of static chaos in spin glasses is studied. For the prob- 
lem of chaos with magnetic field, scaling relations in the case of the SK 
model and short-range models are presented. Our results also suggest 
that if there is de Almeida-Thouless line then it is similar to that of 
mean- field theory for d = 4 and close to the h = axis for d = 3. We 
estimate the lower critical dimension to be in the range 2.7 — 2.9. Nu- 
merical studies at d = 4 show that there is chaos against temperature 
changes and the correlation length diverges like ^ ~ (T — T*)~^. 







The nature of the spin-glass phase is poorly understood in short-range 
spin glasses[jl|. One very interesting topic is the static chaos problem. By this 
we mean how the free energy landscape is modified when a small perturbation 
is applied to the system. 

Most generally one is interested in the problem of chaos when a small 
magnetic field is applied or when the temperature of the system is slightly 
changed. The interest of this problem is twofold. On the one hand, it is 
important to discover what is the nature of the spin-glass phase in short- 
range models. The study of chaos can give interesting predictions regarding 
this question. On the other hand, it is relevant for the understanding of 
some recent cycling temperature experiments in spin glasses [0, ^. Under 
the hypothesis that in the dynamical experiment one is probing some kind of 
equilibrium states |^ it is of the utmost importance to understand the effect 
of changing the temperature on the free energy landscape. 

In this letter we shall present results on the problem of chaos in a magnetic 
field and will see how scaling arguments may be used within the spin-glass 
phase. From these arguments one can predict the nature of the spin-glass 
phase in finite dimensions. Also in the case when the temperature is changed 
we present some results which show that in finite dimensions the system is 
more chaotic to temperature changes than in mean-field theory. 

The study of chaos with a magnetic field was adressed by I. Kondor in 
case of mean- field theory ||^. Let us suppose two copies of the same system 
(i.e. with the same realization of bonds): one at zero magnetic field, the other 
one at finite magnetic field h (in the general case one could study different 
non-zero magnetic fields but for simplicity we will focus on the particular 
case in which one of the magnetic fields is zero.) 

The hamiltonian can be written as: 

i<j i<j i 
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One defines the order parameter 



^ = ^E^^^-^ (2) 



and its corresponding correlation function 



C{i - i) = {(TiajT.Tj) (3) 

wliere (■) means average over disorder and ((■)) is tlie usual tliermal average 
over tlie Hamiltonian eq.(|l|). 

It lias been shown by I. Kondor that this correlation function decays to 
zero at large distances with a finite correlation length which diverges for 
h —>■ 0. For a finite field g = is a stable solution and the system is chaotic. 
In fact, the propagator G{p) (i.e. the Fourier transform of C{i — j)) can be 
exactly computed within the Gaussian approximation. Its singular part is 
given by @]: 

rq,na. rQma. p2 ^ 1 + ^ (g) A ( Q) 



with 



PQmax 

X{q) = - Qmax + / dqx{q)) (5) 

Jq 

where (3 is the inverse of the temperature. The same expression applies in 
the case of \{Q)- Here q{x) and Q{x) are the order parameter functions for 
the spin glass at zero and h field respectively. 

It was found that the correlation length ^ diverges like (1 — \{Qmin))~^ ■ 

2 2 

Close to Tc = 1 we have Qmin ~ hs. This gives ^ ~ /i~3. These are the 
results already obtained in reference 0. 

Now we want to extract more information on the transition when /i — 0. 
We can define a certain kind of non-linear susceptibility by: 

Xni = E (6) 
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It can be shown that it is also given by Xni = N (g^) with q given in eq.(0). 
In terms of eq.(^ we have Xni = G{0). Using the known expressions for 
q{x) and Q{x) close to Tc in eq.(H) we obtain a divergent expression for G{0). 
Its most divergent part is given by 

fQmax rlf) 

G(0) ~ / -^^-^ (7) 

(1 - A(g))t ^ ' 
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which gives Xni ~ ^ ~ h~^. In the Gaussian approximation this gives 
G{x) ~ ^ for h ^ with = d — 4. This means that correlation functions 
decay more slowly in the spin-glass phase than they do at the critical point 
{fi = d — 2) and this is a consequence of the existence of a large number of 
marginal states in the spin-glass phase. 

In general we can introduce exponents 7 and u in the spin-glass phase 
such that Xni ~ h~"' and ^ ~ h~'^. Now we are interested in deriving some 
finite-size scaling relations in order to establish (if it exists) the upper crit- 
ical dimension below which we expect scaling to be satisfied. Also scaling 
relations can be very helpful to test predictions using numerical simulations. 

The derivation of scaling relations within the spin-glass phase proceeds 
analogously as is done at the critical point. At the critical point (below 
six dimensions) we can write Xni ~ L?'~'^ f {N h? q) because h? J2a<b Qat is 
the singular part of the free energy per site {Qat = {o^aCb) with a, b replica 
indices). At the critical point {Qat = q) we introduce the exponent 6 such 
that q = hi . From these results one gets ^ ~ h ^^1^ and the scaling relation 

Xni = L'-V{Lh'^) (8) 

This gives the usual hyperscaling relations (3{5 + 1) = du and S = ^rf^- In 
mean-field theory 0] we have q ^ h ,i.e 5 = 2, z/ = |. Hyperscaling relations 
give rj = and d = 6 which is the upper critical dimension. 

In the spin-glass phase, the derivation is slightly different because Qab 7^ 0. 
To obtain the singular part of the free energy one has to substract from 
h"^ J2a<b Qab that part corresponding to zero magnetic field. In mean-field 
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theory this is given by ^^(/q Q{x)dx — /g'""""' q{x)dx) which is proportional 
to h'^ Qmin Xmin with Xmin squal to the first breakpoint of the function Q{x). 

2 

Because Qmin ~ Xmin ~ hs we obtain 

Xra=NlfiNhf) (9) 
If an upper critical dimension exists then it should be (i„ = 5 because for 

2 

that dimension ^ ~ h~3. 

We have performed Monte Carlo numerical simulations of the SK model 
in order to test the prediction eq.(|^). The results for different magnetic 
fields ranging from 0.2 up to 1.0 at T = 0.6 are shown in figure 1 for several 
sizes (up to = 2016). There is agreement with the prediction eq.(^) and 
Xni versus the field h does not change its behaviour when crossing the de 
Almeida-Thouless line (/iat ~ 0.45 for T = 0.6). 

Now we present the results of our simulations and our predictions for 
short-range models. We expect that scaling is satisfied in the spin-glass 
phase below five dimensions. In general, one has: 

Xnl = fiL^'h^ Qmin Xmin) (10) 

The critical point is a particular case of eq.(|l^). Since there is not replica 
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symmetry breaking one has Xmin ~ 1 and Qmin = q ^ hs . Putting X = 2 — r] 
one recovers eq.(|^). To test this expression at the critical point we have 
performed Monte Carlo numerical simulations of the 4d ±J Ising spin glass 
(with periodic boundary conditions) which is known to have a transition at 
T = 2.06 ± 0.02 @, I. In this case the known values r] ~ -0.25, 6 ~ 3.6 
fit the data reasonably well but even though the scaling is very sensitive to 
the precise value of the critical temperature. This shows that, at the critical 
point, the critical behavior of eq.(|^) is the same as that of the correlation 
function of two identical copies of the system at the same field h. 

Now we present our results at T = 1.5 in the spin-glass phase for the 4d 
Ising spin glass. One has to be sure that samples are well thermalized. To 
this end we have performed a simmulated annealing algorithm which reaches 
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equilibrium in a reasonable time. Our results in figure 2 fit well a scaling law 
Xni ~ L^-^^ f{hL^-^^). This gives A ~ 3.25, ly ~ 0.69 and 7 = Az/ ~ 2.24. We 
should draw attention to the fact the value found for A is close to that found 
in mean-field theory eq.(|^) putting = L^. 

We have also studied the 3d ±J Ising spin glass which has a transition 



close to T = 1.2 |T0|, [TT]. Simulations in the critical temperature give ex- 
ponents in eq.@ in agreement with those already known. We have also 
performed simulations for small sizes at T = 0.8 in the spin-glass phase. Our 
results are compatible with A ~ 2.4, u ~ 0.7 and 7 = Az^ ^ 1.9 

From our results at d = 3, 4 it seems that A = y and = | are a good 
approximation to the exponents at least for d < 5. 

Now we can adress the question of the existence of a phase transition 
line in finite magnetic field in the 4d Ising spin glass (the so called AT line 



[0). From a theoretical point of view the problem remains unsolved [|T5 
Recent numerical studies suggest that this line really exists |T3|, ||, If this 
is the case and there is also mean-field behaviour in the spin-glass phase at 
zero magnetic field |jT6|, l7| then it is natural to suppose that (as happens in 
mean-field theory) Qmin ~ Xmin- This is a very plausible hypothesis which 
agrees with the fact that q{x) ~ x for small x, or that P{0) = (^)q=o has 
a finite value [l^. If Qmin ^ hi (5 = 3 in the mean-field case) we obtain, 

2(2+S) 

from eq.(|l(]|) the result ^ ~ h For d = A this gives 5 ~ 5.27 and 

Qmin ~ hP'^^ . In three dimensions 5 is uncertain but is a very large value (of 
order 20). 

We expect the transition line in magnetic field to occur when Qmin ~ 
Qmax, i-e. /i« ~ or /i ~ with t = Tc — T. Because (3 ~ 0.6,0.5 
in d = 4||l6|,3|T^ respectively this gives h ~ r^'^^ in four dimensions and 
h with r of order five in three dimensions. This means that in three 
dimensions the AT line is very close to the /i = axis and very difficult to 



see numerically, at least not very far from Tc [|T^, |2^ . 

From the scaling relation eq.(lO) we can also estimate the lower critical 
dimension. If we define the exponent 9 as that exponent for which {^'^h'^q) ^ 
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where q is finite then this gives the thermal exponent introduced in droplet 



models ^ The exponent vanishes for d = di where di is the 
lower critical dimension. We have obtained 9 = 1,0.55,0.05 in d = 5,4,3 
respectively. Extrapolating to 5 = we estimate di = 2.7 — 2.9 which is 



in agreement with perturbative calculations in the spin-glass phase ||2^ but 
higher than the value reported in |^ . 

We have also investigated the problem of chaos against temperature 
changes. The outline of the ideas follow that presented above in the case 
of a magnetic field. Now one couples two copies of the system at different 
temperatures. In mean-field theory the problem has not yet been fully solved 
and chaos could be marginal In finite dimensions a interesting be- 

haviour is expected in low dimensions Our numerical results for the SK 
model indicate that, if there is chaos, it is very small (details will be presented 
elsewhere). We have performed simulations in the Ad Ising spin glass. Figure 
4 shows the non-linear susceptibility defined in eq.(^ using eq.(g) which is 
the overlap obtained by coupling two identical copies of the system at differ- 
ent temperatures. Our results are consistent with a correlation length which 
diverges like (T — T*)~^ where T* is the reference temperature of one of the 
two copies. This results are in agreement with perturbative calculations in 
the range of dimensionalities 6 < d < 8 |2^ . 

Summarizing, in the case of chaos with a magnetic field we find that 
there is scaling behaviour in the spin-glass phase in mean-field theory, the 
main result being eq.(|^. Short-range systems also satisfy scaling relations 
from which we can extract the exponents associated to the correlation length. 
We derive that if there is an AT line then in d = 4 it is similar to that of 
mean-field theory and in = 3 it is very close to the h = axis and more 
difficult to see using numerical simulations. The lower critical dimension is 
also predicted to be in the range di ~ 2.7 — 2.9. We also reported some results 
of chaos in temperature which show that short-range models are more chaotic 
than the SK model and the correlation length diverges like {T — T*)^^. 
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FIGURE CAPTION 



Fig. 1 Chaos with magnetic field in the SK model at T = 0.6. Field 
values range from h = 0.2 up to /i = 1.0 for the smaller sizes and up to /i = 0.4 
for the largest ones. The number of samples range from 200 for A'^ = 32 down 
to 25 for = 2016. 

Fig. 2 Chaos with magnetic field in the Ad Ising spin glass at T — 1.5. 
Magnetic field values range from h — 0.1 up to h — 1.. The number of 
samples is approximately 100 for all lattice sizes. 

Fig. 3 Chaos with temperature changes in the Ad Ising spin glass. The 
reference temperature of one copy is T* = 1.5. Temperature values of the 
other copy range from 1.6 up to 2.2. The number of samples is approximately 
100 for all lattice sizes. 
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